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ABSTRACT 

We describe the discovery of HAT-P-4b, a low-density extrasolar planet transiting BD+36 2593, 
a V = 11.2 mag slightly evolved metal-rich late F star. The planet's orbital period is 3.056536 ± 
0.000057d with a mid-transit epoch of 2, 454, 245.8154 ± 0.0003 (HJD). Based on high-precision pho- 
tometric and spectroscopic data, and by using transit light curve modeling, spectrum analysis and evo- 
lutionary models, we derive the following planet parameters: M p = 0.68±0.04Mj, R p = 1.27 ± 0.05 Rj, 
p p = 0.41 ±0.06 g cm" 3 and a — 0.0446±0.0012 AU. Because of its relatively large radius, together with 
its assumed high metallicity of that of its parent star, this planet adds to the theoretical challenges 
to explain inflated extrasolar planets. 

Subject headings: planetary systems: individual: HAT-P-4b — stars: individual: BD+36 2593 - 
techniques: photometric — techniques: spectroscopic 



1. INTRODUCTION 

In the course of our ongoi ng wide field plane tary tran- 
sit search program HATNet (jBakos et al.ll2004f ). we have 
discovered a large radius and low density planet orbit- 
ing an 11th magnitude star BD+36 2593. This planet 
is the fifth member of a group of low-density transiting 
exoplanets. The combination of its low mass and the rel- 
atively high metallicity and age of the parent star makes 
theoretical interpretation of its large radius difficult. In 
this Letter we describe the observational properties of 
the system and derive the physical parameters both for 
the host star and for the planet. We also briefly comment 
on the theoretical status of inflated extrasolar planets. 

2. THE PHOTOMETRIC DISCOVERY AND FOLLOW-UP 
OBSERVATIONS 

The star BD+36 2593 (also GSC 02569-01599) at 
a = 15' 1 19 m 57 ! r92, S = +36°13'46'/7, is contained in field 
G191 of HATNet, centered at a = 15^28™, 5 = +37°30'. 
As we show in the remainder of the paper, the star is 
orbited by a planetary companion, and so we label the 
host star as HAT-P-4 and the planet as HAT-P-4b. Field 
G191 was monitored from 2004 December until 2005 
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Fig. 1. — Folded and unbinned light curves of HAT-P-4. The 
HATNet and FLWO 1.2 m data are plotted in the upper and lower 
panels, respectively. The transit model fit to the FLWO data 
is shown by continuous line (see Table [3] for the resulting transit 
parameters). The out-of-transit level of the HATNet data corre- 
sponds to I = 10.537 mag, phase is zero at the transit center. 

June by the HAT-7 telescope at the Fred Lawrence Whip- 
ple Observatory (FLWO) of the Smithsonian Astrophys- 
ical Observatory (SAO). Nearly 90% of the data points 
were gathered after February 2005. The field is relatively 
sparse, with ~ 14000 objects brighter than / rj 14 mag. 
Most of the light curves have some 5300 data points. 
The sampling cadence is 5.5 min, slightly longer than the 
5 min integration time. 

A first look at the compact subset collected in the 
98 d time span between March and June, 2005, re- 
vealed the presence of a transit signal with a nearly in- 
teger 3 d period. The detection was made possible by 
the application of the Trend Filtering Algorithm (TFA; 
iKovacs. Bakos. fc Novesll2005[ ). However, assuming the 
primary is a main-sequence F or G type star, the rela- 
tive length of the transit is some 20 — 30% greater than 
expected for a Jupiter-size companion with the above or- 
bital period. Although both the lack of detection in the 
raw data and the length of the transit made us suspicious 
about the viability of this candidate, we left some room 
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for the possibility that the primary was slightly evolved 
(an assumption that proved to be true in the subsequent 
investigations — see §S]). Therefore, we proceeded with 
rejection-mode spectroscopy, which showed no sign of ra- 
dial velocity (RV) variation at the kins" 1 level. 

An improved reduction of all available frames of the 
field was completed by February 2007. In addition to 
using a new data pipelin e that employs refined astrom- 
etry (|Pal fc Bakosl [20061 ) and aperture photometry, we 
detrended the light curves before signal search with the 
aid of an External Parameter Decorrelation technique 
(EPD, see also Bakos et al. 2007). This technique uti- 
lizes the fact that various "external parameters" that are 
specific to the star, such as sub-pixel position on the 
frame, point spread function properties (e.g., width and 
elongation), or specific to the frame, such as telescope 
position, are correlated with the deviations of the star's 
brightness from the median. The technique is optimal 
for stars that are otherwise non- variable most of the time 
(e.g. transit candidates). EPD derives the correlation be- 
tween brightness deviations and the underlying external 
parameters, and subsequently corrects them for each in- 
dividual star. This is different from what is done during 
the application of TEA, where we consider the full time 
history of the light variation and use the light curves to 
"cure themselves" by recognizing the hidden systematics 
in each other. The two methods are complementary and 
we use both (EPD followed by TFA). All these led to a 
powerful confirmation of the earlier detection. 

The folded light curve constructed from the HAT- 
Net data is shown in the upper panel of Fig. [TJ 
Th e best period was obtained fr om the BLS analy- 
sis ([Kovacs. Zucker. fc Mazehll2002l ) after applying TFA 
with some 700 template stars. 11 As seen, the transit is 
clearly visible in the light curve even though no binning 
was applied. There is no sign of periodic out-of-transit 
variation, indicating that this is a near "textbook" tran- 
sit signal. 

To update the photometric ephemeris and obtain a pre- 
cise light curve for model fitting, we observed our target 
with KeplerCam on the 1.2 m telescope at FLWO. We 
observed two full transits on May 24 and 27, 2007. The 
folded light curve in the Sloan z band for all 985 data 
points is shown in the lower panel of Fig. [T] The follow- 
up observations verify the discovery light curve and are 
accurate enough to make transit modeling possible. 

3. THE SPECTROSCOPIC VERIFICATION AND 
EXCLUSION OF BLEND SCENARIOS 

As mentioned in § [2j basic rejection mode spectroscopy 
to exclude stellar binary and certain blend types was 
already initiated after the first detection in the fall of 
2005. By June 2006 we had collec ted 9 spectra with 
the CfA Digital Speedometer (DS; iLathaml I1992D on 
the FLWO 1.5 m telescope, covering a 45 A window cen- 
tered at the Mg b triplet at 5187 A. We derived spec- 
troscopic parameters by comparing the observed spec- 
tra against synthet ic spectra, as described in detail by 
iTorres et~a l. (2002). The initial stellar parameters from 
this analysis were: T e g= 5500 K, log 5= 3.5 [cgs] and 
vsini= 5.9kms , indicating an evolved G-type pri- 

11 The EPD technique played the main role in the detection; 
TFA subsequently led to an increase of 20% in the signal-to-noise 
ratio. 



TABLE 1 

HIRES Relative Radial Velocities for HAT-P-4. 



BJD 


RV 


CTRV 


O-C 


Phase a 


(2,400,000+) 


(ms ) 


(m s J 


(m s ) 




54186.98522 


101.0 


2.1 


2.7 


0.753 


54187.11241 


92.6 


2.1 


-2.6 


0.794 


54188.01160 


-25.4 


2.0 


0.3 


0.088 


54188.07150 


-35.8 


2.0 


-2.0 


0.108 


54189.00174 


-28.1 


2.2 


-2.8 


0.412 


54189.08262 


-13.9 


2.0 


-0.7 


0.439 


54189.13222 


-3.1 


3.3 


2.3 


0.455 


54249.93769 


-52.5 


2.8 


-3.6 


0.349 


54279.86357 


-36.8 


2.9 


8.4 


0.139 


a Relative to the 


center of t 


ransit. 








0.4 0.6 
Photometric phase 

Fig. 2. — Relative radial velocity and bisector span variations 
from the Keck observations of HAT-P-4. Panel (a) shows the nine 
Keck RV measurements with a zero-eccentricity orbital fit. In panel 
(b) we show the bisector spans (from nine iodine exposures and one 
template spectrum) displaying a much smaller scatter. Error bars 
include the 3 ms — 1 estimated velocity jitter. 

mary 12 . The DS spectra were also used to derive ra- 
dial velocities by cross-correlating the spectra with syn- 
thetic sp ectra based on Kurucz model atmospheres (e.g. 
ILathaml [2002). The rms scatter of the RV data was 
0.71 kms _1 around the mean velocity of — 1.58 kms -1 . 
The velocities did not show any correlation with the 
orbital phase and we found no sign of a second stel- 
lar component in the spectra. These spectroscopic re- 
sults, together with the secure detection in the HATNct 
data, made a viable case for obtaining high-precision RV 
measurements to seek evidence for orbital motion and 
to obtain refined stellar parameters. Accordingly, high- 
resolution s pectroscopy was conducted with the HIRES 
instrument (jVogt et aLlll994h on the Keck I telescope be- 
tween 2007 March 27 and May 29. See, e.g.. lTorres et all 
(|2007f ) for the details of the observational procedure. The 
nine resulting RV measurements are listed in Table [TJ 
They are relative velocities in the Solar System barycen- 
tric frame of reference (see iButler et al.[|1996h . 

12 These parameters were derived with [Fe/H]= 0.0. By using 
[Fe/H]= 0.24 from the Keck spectra we get a slightly better match 
with the synthetic spectra that yields T c ff= 6000 K and logg= 4.0, 
in better agreement with the Keck results described below. 
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The data were fitted by a circular 13 orbit with period 
and mid transit epoch constrained by the photometric 
data (see Table [3]). The result of this two-parameter 
(systematic offset and semi-amplitude) fit is displayed 
in Fig. [21 The unbiased estimate of the standard devia- 
tion of the residuals is 4.1ms _1 , larger than the internal 
errors listed in Table [TJ A part of this scatter may come 
from the last data point that was taken under unfavor- 
able weather conditions (leaving out this point we get 
2.7ms" 1 for the standard deviation of the residuals). 
Nevertheless, it is also possible that there is a "veloc- 
ity jitter" due to stellar activity. The size of this jitter 
is estimated to be about 3ms" 1 , fall ing in the exp ected 
range for an inactive late F-type star (Wright 2005) . The 
low activity level is supported by the absence of emission 
features in the Ca II H and K lines. 

A crucial part of the verification of the sub-stellar na- 
ture of the companion is testing various blend scenar- 
ios. There are at least two basic ways of such testing: 

(i) light curve modeling, combined with the information 
available from the spectra (e.g., maximum flux contri- 
bution by the blended binary); (ii) checking spectral line 
asymmetry by searching for bisector sp an variations (e.g., 
ISantos et all 120021 : [Torres et al.ll2005D . Our experience 
shows that while (i) requires several pieces of informa- 
tion (good quality full light curve, running many binary 
models and performing a fairly deep spectrum analysis), 

(ii) is more sensitive to hidden components and simpler 
to perform. Therefore, we settled on (ii) and searched 
for a variation in the differences measured between the 
velocities at the top and the bottom of the correlation 
profiles of the individual spectra. The result is shown 
in the lower panel of Fig. O For a hidden stellar binary 
blended by our target we would expect a variation in 
the bisector span in phase with the radial velocity, and 
with an amplitude c omparable to that of the measured 
RV variation itself ([Torres et al.l l2005h . However, the 
standard deviation of the bisector variation is 6.5 ms" 1 , 
which is only 8% of the RV amplitude (the same num- 
bers are 2.9ms" 1 and 4%, respectively, if we leave out 
the last data point). Furthermore, no correlation with 
the RV variation is seen. Therefore, we are confident 
that the source of the RV variation is stellar wobble due 
to the gravitational pull of a sub-stellar companion. 

4. STELLAR AND PLANETARY PARAMETERS 

The cornerstone of the derivation of the absolute pa- 
rameters of a planet discovered by radial velocity and 
transit observations is the accurate estimation of the stel- 
lar mass and radius. 

The procedure often involves determination of T c g and 
[Fe/H], as well as logg, from high-resolution spectro- 
scopic analysis. If an accurate distance is available, for 
example from a Hipparcos parallax, then the absolute 
magnitude can be used to improve the value for the stel- 
lar radius. In general, none of the methods yield bet- 
ter than ~ 10% accuracy in the derived stellar param- 
eters (implying corresponding limits to the accuracy of 
the absolute planet parameters). An improvement can 
be achieved by using the value of a/ J?* derived from the 
transit light curve to measure directly the mean stellar 



TABLE 2 
Stellar parameters for HAT-P-4. 



Parameter 



Value 



Source 



T cS (K) 
[Fe/H] (dex) 
v sini (kms ) 
Mass (M©) 
Radius (Rq) 
logg (cgs) 

L* (£©) 
My (mag) 
Age (Gyr) 
Distance (pc) 



5860 ± 80 
+0.24 ± 0.08 
5.5 ± 0.5 



1.26 
1.59 
4.14 
2.68 
3.74 
4.2 [ 
310 ± 30 



-0.14,0.06] 
-0.07,0.07] 
-0.04,0.01] 
-0.34,0.39] 
-0.16,0.16] 
-0.6,2.6] 



SME a 

SME 

SME 

Y 2 +LC+SME b 

Y 2 +LC+SME 

Y 2 +LC+SME 

Y 2 +LC+SME 

Y 2 +LC+SME 

Y 2 +LC+SME 

from M v c 



a SME — "Spectroscopy Made Easy" package to gen- 
erate synthetic spectra and to fit the observed ones 
iValenti fc Piskunov|ll996l) , b Y 2 +L C+SME = Yonsei-Yale 
isochroncs (Dcmarqucct al. 2004). transit light curve 
modeling and SME analysis. Using V = 11.22 ± 0.12 of 
IDroegel J2006D and assuming zero reddening. 



13 A fit with free-floating eccentricity yields e 
supporting the assumption of a circular orbit. 



0.055 ± 0.035, 



density (|Sozzetti et alj |2007) . Since this latter parameter 
is directly linked to log g, entering in the spectroscopic 
analysis, we followed an iterative determination of the 
stellar and planetary parameters. 

First, for the determination of [Fe/H] and T c g, we 
used the iodine-free template spectrum from Keck. 
The modeling was performed using the SME soft- 
ware ([Valenti fe Piskunovl 1 199(f ) incor porating the same 
metho d and atomic data as given in IValenti fc Fischerl 
(2005). We obtained T off = 6032±80K, [Fe/H]= +0.32± 
0.08 [dex] and logg= 4.36 ± 0.11. 

Next, for the computation a/i?*, we fitted the high 
precision K eplerCam light curve (s ee Fig.[T]) by using the 
formulae of iMandel fc AgoU (120021) with quadratic limb 
darkening coefficients from lciaretl ( 2004D . We set e = 0.0 
as a result of our test on the RV data with non-zero 
eccentricity. Fitted parameters were the center of transit 
T c , the radius ratio Rp/R*, the normalized relative semi- 
major axis a/ R±, and the impact parameter b. 

Next, to compute the stellar mass and radius we re- 
lied on current stellar evolution models. As in our ear- 
lier papers, we compared the observational properties of 
the host star w ith a finely inter p olated grid of model 
isochrones from iDemarque et al.l (|2004[ ). The inferred 
stellar properties are based on the best match to the 
measured values of T e ff, [Fe/H], and a/R+ within their 
observational errors, in a \ 2 sense. This procedure led 
to a better approximation of the stellar gravitational ac- 
celeration with log 3= 4.12 ± 0.04. 

In the second loop of iteration we used the newly deter- 
mined value of log g in the SME analysis and redid the 
above sequence of computation. This led to a slightly 
modified set of stellar and planetary parameters with 
somewhat smaller errors than in the first loop. Since 
the changes were, in general, fairly small (e.g., logg has 
changed to 4.14± 0.03), we decided to stop the iteration 
after this second loop. The final stellar and planetary pa- 
rameters are shown in Table [2] and Table [3j respectively. 

Concerning the derived parameters and their errors, 
we note the following. The dependence of the result 
on the evolutionary models was tested by using the 
isochrones given for solar-scaled Z = 0-03 m odels with 
core overshooting bv lPietrinferni et~a l. (2006). As noted 
by the authors, their models are hotter by some 200 K 
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TABLE 3 

Spectroscopic, light curve and planet 
parameters of hat-p-4 



P a T*ti T'ri ptpr 


Value 


Spectroscopic parameters: 




Period (d) a 


3.056536 ± 0.000057 


T c (HJD) a 


2,454,245.8154 ± 0.0003 




2,454,248.8716 ± 0.0006 


K (ms" 1 ) 


81.1 ± 1.9 


Offset velocity (ms -1 ) . 


12.1 ± 0.9 


e c 


0.0 


Light curve parameters: 




Transit duration (day) . . . 


0.1760 ± 0.0003 


a/ 


6.04 [-0.18,0.03] 


Rp/FU 


0.08200 ± 0.00044 


b = a cos i 1 R* 


0.01 [-0.01,0.23] 


Planet parameters: 




M p (Mj) 


0.68 ± 0.04 


Rp(Rj) 


1.27 ±0.05 


Pp (gem -3 ) 


0.41 ± 0.06 


logg p (cm s -2 ) 


3.02 ± 0.02 


a (AU) 


0.0446 ± 0.0012 


V(deg) 


89.9° [-2.2,0.1] 



a Taken from the photometry (HATNet and 2 nights of 
FLWO 1.2 m ). b The 7 velocity is -1.58 ± 0.24kms _1 , 
from the DS data. c Adopted 



than those of iDemarque et al.1 (|2004f ). Therefore, we 
used an effective temperature of 6060 K. With these 
input parameters we got nearly the same stellar and 
planet parameters as from the Yonsei-Yale isochrones 
(i.e., p p = 0.41 gem -3 , age= 4.5 Gyr). By using models 
without overshooting or of lower effective temperature, 
we got larger ages and also slightly larger densities, up 
to p p = 0.43 gem -3 . The stability of the planet den- 
sity is mainly related to the correlated change of stellar 
mass and radius when models or input parameters are 
changed. These age ranges and the derived metallicity 
fit rea sonably to the relation recently given bv lReid et al.l 
(|2007D . We also note that the derived radius is 34% 
larger than the one corresponding to an unevolved main- 
sequence star. This explains the longer than expected 
transit duration by which we were puzzled at the early 
phase of the discovery. 

5. DISCUSSION AND CONCLUSIONS 

We presented the discovery data and derived the phys- 
ical parameters of HAT-P-4b, an inflated planet orbiting 
BD+36 2593. Among the 20 transiting planets discov- 
ered so far, there are five with p p < 0.4 g cm -3 . All others 
have at least 50% higher densities. For ease of compar- 
ison, Table 0] lists the relevant properties of the five in- 
flated planets. It is remarkable how similar these planets 
are (except for TrES-4 that has distinctively low density) . 
With its Safronov number of 0.036, HAT-P-4b belongs to 
th e Class II planets accordin g to the recent classification 
of lHansen fc Barman! (|2007f ) and (together with TrES- 
4) further strengthens the mysterious dichotomy of the 
known transiting planets in this parameter. The par- 
ent star of HAT-P-4b is among the largest radii, largest 
mass, lowest gravity and highest metallicity transiting 
planet host stars. 



Current models of irradiated giant planets are able to 
match the observed radii of most of the planets with- 



TABLE 4 

Comparison of the properties of inflated planets. a 



Name 


P 


a 


M 


R 


P 


logs 




(d) 


(AU) 


(Mj) 


(Rj) 


(CGS) 


(CGS) 


WASP- lb 


2.52 


0.038 


0.87 


1.40 


0.39 


3.04 


HAT-P-4b 


3.06 


0.045 


0.68 


1.27 


0.41 


3.02 


HD 209458b 


3.53 


0.045 


0.64 


1.32 


0.35 


2.96 


TrES-4 


3.55 


0.049 


0.84 


1.67 


0.22 


2.87 


HAT-P-lb 


4.47 


0.055 


0.53 


1.20 


0.38 


2.96 



a Data from IShporer et all (|2007). this paper, Burrows et al.l 
120071) . iMandushev et alj (t200W l and lWinn et al.l J2007D . From 
top to bottom, metallicities for the parent stars are: 0.23 
IStempels et al.N2007D . 0.24, 0.02 0.0 (adopted) and 0.13. 



out invoking any additional heating mechanism. Higher 
metallicity cases, such as the present one, however, may 
pose problems (assuming that the planet and star have 
similar metallicities) . More metals imply two opposite ef- 
fects on the radius: (i) inflating it due to higher opacities 
in the envelope; (ii) shrinking it due to the higher molec- 
ular weight of the interior and the possible development 
of a large high den sity core. These effect s have been dis- 
cussed recently by iBurrows et al.1 (|2007f ). Since WASP- 
lb is similar in several aspects (i.e., irradiance, metal- 
licity) to HAT-P-4b, we consid er the coreless model s of 
WASP-lb as shown in Fig. 7 of lBurrows et all (|2007f) . It 
seems that HAT-P-4b can be fitted by near solar metal- 
licity coreless models, assuming that its age is not too 
much greater that 4 Gy r. We also refer to the la yered 
convective mechanism of lChabrier fc Baraffd (|2007f) that 
gives an alternative explanation for planets with inflated 
radii. 

We conclude that more definite statements on the re- 
lation of the observations and planet structure theories 
can be made only by reaching higher accuracy in the ob- 
served star/planet parameters. Nevertheless, HAT-P-4b 
(together with WASP-lb) does not seem to support the 
existence of a simple relation betw een host star metal- 
licity and planet's core mass (see iGuillot et all 120061 : 
IBurrows et aHl2007t ). 
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